Guchi-TP-003 



Exact Solutions for Boson- Fermion Stars in (2+1) dimensions 

Kenji Sakamoto^ and Kiyoshi Shiraishi^'^* 
^ Graduate School of Science and Engineering, Yamaguchi University 
Yoshida, Yamaguchi- shi, Yamaguchi 753-8512, Japan 

"^Faculty of Science, Yamaguchi University 
Yoshida, Yamaguchi- shi, Yamaguchi 753-8512, Japan 
(February 3, 2008) 



Abstract 

We solve Einstein equations coupled to a complex scalar field with infinitely 
large self-interaction, degenerate fermions, and a negative cosmological con- 
stant in (2 + 1) dimensions. Exact solutions for static boson-fermion stars 
are found when circular symmetry is assumed. We find that the minimum 
binding energy of boson-fermion star takes a negative value if the value of the 

cosmological constant is sufficiently small. 
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I. INTRODUCTION 



Lower-dimensional gravity provides us with a useful arena where we can often construct 
exactly solvable models and can study the model analytically. The solutions for relativistic 
boson stars ( for reviews) are only numerically obtained in four dimensions. In the 
previous paper ||^, we obtained an exact solution for a nonrotating boson star in (2 + 1) 
dimensional gravity with a negative cosmological constant. We assumed that the scalar field 
has a strong self-interaction. An infinitely large self-interaction term in the model leads to 
much simplication as in the (3 + 1) dimensional case [Q. 

In the present paper, we derive an exact solution for a nonrotating boson-fermion star 
in (2 + 1) dimensional gravity with a negative cosmological constant. We assume that the 
scalar field has a strong self-interaction as in the previous paper. 

Four-dimensional boson-fermion stars have been studied by many authors . Our three- 
dimensional model includes an additional constant, the cosmological constant. We will see 
that the value of the cosmological constant is restricted by the stability condition of the 
boson-fermion star. We will also find, nevertheless, that some physical quantities of the 
boson-fermion star behave as in the four-dimensional case. 

In Section ^ we show the field equations which we solve. In Section |T|, we obtain exact 
solutions for the field equations shown in Section We pay our attention on the expression 
of mass of the boson-fermion star. The particle numbers of bosons and fermions for the 
boson-fermion star described by the exact solutions is given in Section In Section 0, we 
restrict our attention to the fermion star, i.e., the case that the particle number of bosons 
vanishes. We discuss stablity of the boson-fermion star in Section^. Section [VI I| is devoted 
to conclusion. 
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II. FIELD EQUATIONS 



We consider a complex scalar field if with mass ttlb and a quartic self-coupling constant 
A. We assume that the coupling constant takes an infinitely large value. The equation of 
motion for the scalar field is 

V^V9 — m|y9 — A Iv^l^ (y9 = 0. (1) 

In addition, we assume that a degenerate fermion field with mass mp is coupled to 
gravity. In this paper, we assume no direct coupling between the boson and fermion fields. 
The Einstein equation is written as 

K - = SvrG (t^I: + T^^) + C6i:, (2) 

where the positive constant C stands for the (negative) cosmological constant. G is the 
Newton constant. 

The energy-momentum tensor of the boson field is 

T^^; = 2Re ( W V.^ - I \Vvf 5^;) - m| 6^^ -^\^\' 5^ (3) 

The energy-momentum tensor of the fermion field is written in the perfect fluid form as 

T^^^ = dzag.{-p,P,P), (4) 

where the hydrostatic equilibrium is assumed. 

In three dimensions, the energy density and pressure are given by 

and 

respectively. Here /i stands for a chemical potential for the fermion field. 
The number density of fermions is given by 



We assume the three-dimensional metric for a static circularly symmetric spacetime as 



-a\r)dr + (3\r)dr' + 



2 , 



d9' 



(8) 



where a, (3, and 7 are functions of the radial coordinate r only. The constant 70 will be 



determined later (see Eqs. (|35| - p6|) and below.). Because of the coordinate invariance with 
respect to r, we have a freedom in choosing the function /5(r). We will use this residual 
"gauge choice" later. 

We also assume that the complex scalar has a phase which is linear in the temporal 
coordinate: 



(9) 



where ^{r) is a function of the radial coordinate r only and is a constant. 

For a static circularly symmetric spacetime, p and P have only dependence on r. Thus 
/i is also assumed to be a function of the radial coordinate r only. 

Then the Einstein equation (0) implies the following equations: 



1 d fld-^\ 
P'y dr \(3 dr ) 
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The equation of motion for the scalar field (0) becomes 

1 d f a^y dip^ 



aP'j dr \ 13 dr 
The equilibrium condition is reduced to 



P 
P 

P 



+ C, 



c. 



0. 



(10) 

(11) 

(12) 



(13) 



To simplify the equations, we make use of rescaled variables as 

A 



f = m^r, (p = VSnGip, A = „ „ y , Cj = u/niB, 
C = C/ml, p = {8nG/ml)p, and P = {8nG/ml)P. (15) 

We also rescale as 

/X = ii/mp. (16) 

In order to study the limit of large self-inteaction, we rescale the variables once more. 
New variables are: 

= f/ a/a, = v^f^, = AC, 

= Ap, and = AP. (17) 
In the limit of large self-coupling, A — * oo. The equations (p^p^) will be reduced to 

i^^^4.:-.;-i.t+p..a, (19) 

>p.x(^-l-vl]^ 0, (21) 



{p^ - l) X [p' + ^p] = 0, (22) 



where ' denotes the derivative with respect to r^,. Here we have discarded the terms of 
order 1/A and higher. Then the terms including ^9^,' disappear in the equations. This 
scheme to study the large self-coupling case was originally adopted by Colpi, Shapiro and 
Wasserman 01. 



III. EXACT SOLUTIONS 



Now, let us solve the set of an algebraic equation and differential equations ( |T^ - |2^) 
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Eq. pip can be solved easily. One solution for is 



^l = -,-l. (23) 

We assume that this solution describes the configuration of the boson field for < r^:B- 
Another solution for (p^, is 

= 0. (24) 

This solution is taken as the configuration of the boson field for r^, > t^b- Note that T^'^ = 
for > r^B- At the boundary r^, = r^B, Oi takes the value a = ui. 
Eq. (|2^ ) can be solved easily, too. One solution for jl^ is 

= 4' (25) 

where aj? is a constant. We assume that this solution describes the configuration of the 
fermion field for r^, < r^:p. Another solution for /i^ is 

fi^ = 1. (26) 

This solution is taken as the configuration of the fermion field for r^, > r^,^. Note that 
T^^^ = 0, or p^, = P^, = for > r^^?. At the boundary = r^^?, a takes the value a = ap. 
From Eqs. (|T^) and (PD|), we find 



a' 



0. (27) 



To solve this, we use the residual gauge choice on /?. For simplicity, we take 

a/37 = Kr,, (28) 

where A' is a constant. Actually, in order to obtain the solution, one can take an arbitrary 
choice. This is because, as we will see later, we will solve the equations by eliminating the 
radial coordinate variable. Then Eq. (^) is reduced to 



where A is a constant. 

From Eqs. and (|28|), we findQ 



X2 r 

Here using Eq. (PW), we obtain 



I W 1 . 
A^ = -<ft + P. + C,. (30) 



Before exploring the exphcit solution for the metric, we rewrite the radial line element 



as 



Consequently, we get 

P'dr' = A#I7^ ^ V2dl'. (33) 

First, we solve the exterior solution for r^, > r^,o, where r^,o = Max[r^:B,r^:F]- Since the 
total energy momentum tensor vanishes for r^: > r^oi the solution of Eq. (|3l]) is found to be 

A^ 



(^7^ - 8GM,) , (34) 



where is an integration constant. 

The full line element for the exterior solution is obtained, by substituting Eqs. ( PB| ) and 
(H) into (D, as 

where a new radial coordinate R is defined as 



^This equation ( |30| ) is consistent with the other equations. 
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We should remember here that C* = Ckjrr? . We take 70 = here to avoid an appear- 
ance of a conical singularity . This choice makes R a usual radial coordinate when Q varies 
in the range < < 27r. 

After rescaling -j^dt dt, we find the metric is precisely the same as the well-known 
BTZ vacuum solution ^ . Thus we identify the constant Mo with the BTZ mass of the star. 

Next we must turn to solving the interior solution of the star. We first consider the 
innermost region, < r^i, where r^i = Min[r^:B,r^:p]. From Eqs. (PB[), and (PH), we 
obtain a differential equation: 

where = (2G'Am^)/(3m|). 

The solution of this equation is given by 

iv V 6*/ 2a^ uj^ 2 \ a J 

(38) 

where is a constant. 

For the intermediate region r^:i < r < r^:o, we must consider two cases separately: for a) 
> r^:p and b) r^:F > r^B- 
Case a): r^B > r^p 

In the region r^:p < r* < r*^, one can solve the equation (^Tj) with (^) and and get 

where Z)a is an integration constant. 

At r=K = r^:p, when a = cti;', the solutions (pHI) and (|39|) must be connected, thus we find 



Da = B. (40) 

The value of B is determined from the condition at the outermost boundary of the star, 
r* = r^^B) where a = Co. From the condition for a smooth connection of the solutions 
and (|39|) , one can find 
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D„ = B = sGMo. (41) 

Case b): r^B > r^p 

In the region r^p < t* < r*^, one can solve the equation (^Tj) with (|^) and (125|) and get 



^ (^y - ^) = ^* (-2^ - 6«F« + + 6«^j + (42) 

where Df, is an integration constant. 

At = r^,B, where a = cD, the solutions ( pS]) and ( 142|) must be connected, thus we find 



= 5. (43) 

The value of B is determined from the condition at the outermost boundary of the star, 
= r=Ki?, where a = ap- From the condition for a smooth connection of the solutions 
and one can find 



( 2 8GMo\ IcD^ 2, «^ 1 2 ^ / "f 2 2\ ^ 2 

' 7 TT^ = -o ^ - In ^ + T^a + i^* -2— - 6aF« + + + ^al 



Db = B = 8GMo. (44) 

In both cases (a and b), we have found that the innermost solution is given by 

^ / 2 8GMo\ IcD^ 2, «^ 1 
i^n^ -^J=-2^-" ^"^ + 2' V « 

(45) 

In the region r=Kj < r,,, < r=Ko, the solutions are given by (|39|) and (^21) with Da = Db = 8GMo 
for Case a) and b), respectively. 

In our solutions, the origin is located at i? = 0, or 7 = 0. At the center of the star, we 
choose the following condition: 

^[^^to + P*o + a) =1, (46) 

where = "(7 = 0), (flo = <^^(7 = 0) = tl^V^o - 1. and P*o = P*il = 0) = F^{a%/al - 
Sap/ao + 2). 

With this choice, one can find that the limit of "no matter" fields yields the three 
dimensional anti-de Sitter space with no conical singularity [||. 

9 



Consequently, we obtain the relation between the BTZ mass of the boson-fermion star 

and v^^o and /iq = ap/aQ from Eq. (^): 

^ \ip% + ip% - (y^,^o + 1) In {ip% + 1) + 2F, (/Ip - 1)^ - q ^^^^ 
(i<^to + ^*(/ig-3/io + 2) + C,)' 

The absence of boson and fermion fields is achieved in the limit ip^o = and jlo = 1- If 
one define the intrinsic mass Mj as 

M, = Mo + ^, (48) 
the intrinsic mass Mi vanishes when in the limit if^o = and /io = 1. 



IV. PARTICLE NUMBERS 

The particle number A^^ of the boson field is given by 

Nb = J (fx^g \g''\i {^dt^* - ip*dtv) ■ (49) 

Therefore we can calculate the particle number of the boson by making use of the above 
solution. It is given by 

/• 2mBU] 2 1 a/57 , 
J SttGA 7o 



2GA70 



J V mi 



B 

- 1 1 —'—da' 



4GA70 V rn% V"^ ) A 

niB I A Kao uj I uP' uP''^ 
^ — 1 — m ^ 



The particle number of fermions is given by 

Up 



Np = d X\/ —g 



V-9tt 



2.1^(^-1] ^dr 
47r \ / 7o 
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(50) 




rripA 



1 3F, {flo - if 



(51) 



One can find that for small (p^Q and jlo — 1, the mass and the partcle numbers are given 
approximately by 



4 Q TP 

Mi ^ + (/2o - 1)' « rriBNB + mpHp. (52) 



V. FERMION STARS 



For a while, we will consider the fermion star, which can be obtained by setting 93*0 = 0. 
In this case, both maxima of the mass and fermion number as functions of </7*o are located 
at /Jo = /iom- P^om is given as 



1 = 



2C* 



1/3 



3C 



1/3 



(53) 



The maximum value of Mj, which we will denote as Mj^ is found to be 



24G 



1/3 



2 ^SG' 



(54) 



The maximum value of Np, which we will denote as Nprn is found to be 



AT, 



The binding energy Ei 



bm 



Fm 



Mi, 



y 3C 



1/3 



(55) 



mpNp^ of the fermion star with maximum fermion 



number is therefore given by 
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If the value of the cosmological constant is larger than Ccru, this state is energetically 
unfavorable because Ehm is positive. The critical value is 

Ccr^t = VSGml. (57) 

It is worth noting that the value of maximum BTZ mass {Mom = — g^) is positive 
for all finite values of C. In particular, in the limit C = 0, the value of Mom vanishes. The 
behavior of Mom for small C is given by Mom oc The value of Npm in the limit C = 

is l/iAGmp). 

If the minimum binding energy is positive, the fermion star is no longer energetically 
favorable. Thus when C > ^/SGm^p, there is no static configuration of the fermion star. 

VI. STABLITY OF BOSON-FERMION STARS 

The stability of four-dimensional boson-fermion stars has been studied by several 
authors. The perturbative stability can be judged by the line of zero eigenvalue about 
particle numbers in the {x, y) plane, where we denote ip^o and jiQ — 1 as x and y, respectively. 
The condition of zero eigenvalue is 

= 0. (58) 

In our solutions, the condition can be reduced to 

x^ (-6^ + Ax^ + x^ + ^F^y^^ + (2C^ - Ax"^ - 3x^ - F^y^^ In (l + x^) = 0, (59) 

where x = ip^:0 and y = fiQ — 1. 

In FIG. [T| we show the line of zero eigenvalue and contour lines of constant Nb and Np 
in the (x, y) plane. The perturbatively stable solutions for boson-fermion stars lie on the 
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dx dy 

dNp dNp 
dx dy 



left of the line of zero eigenvalue. The solutions on the line is absolutely stable, in the sense 
that the stars does not exhibit radial oscillations. 

In our three-dimensional model, the binding energy of the star can take positive value 
even in the left region of the line of zero eigenvalue. Thus the minimum of the binding 
energy is not always a global minimum. 

The contour lines of the binding energy Eb = Mi — tueNb — mpNp is also plotted in 
FIG. |1]. 

The minimum binding energy point is given by the crossing point of zero-eigenvalue line 
(^) and the following line: 



C,y [-2C, + 2x^ + x^ + F,y [6 + 3y + y^ 
^x^ + F,y^ (3 + + a 



+ 3x^ + F^y^ -2{y + 2) xVl + x^ 



(60) 



where x = ip^Q and y = flo — 1. 

Unfortunately, the values of x = ip^:0 and y = jiQ — l which yield the configuration of the 
minimum binding energy cannot be obtained analytically. 

The value of minimum binding energy E^m is plotted as a function of F^ and in FIG. ^. 
One can see that the minimum binding energy becomes positive for sufficiently large and 
small F^. 

The critical line, i.e., the line of zero minimum binding energy is shown in FIG. |^. In 
the limit of F^ = 0, the line ends with C* = C=k crit, where C^, cru ~ 2.5268 is the critical 
value obtained in for boson stars. This fact is trivial because the limit = means no 
contribution of fermions. 

The broken line in FIG. ^ shows the line of 



(61) 



It is obvious that the critical value (|57|) can be found for fermion stars, when the contribution 
of fermions is dominant. 
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If the minimum binding energy is positive, a boson-fermion star with any particle num- 
bers has positive binding energy. Therefore, for arbitrary values of F*, no stable boson- 
fermion star exists for sufficiently large value of C*. The critical value is given approximately 

by 

,rit{F.) ^ 2.5268 + -^F,. (62) 

Finally, we give the mass and particle numbers in a special limit, C* << 1 and F* << 1. 
In this case, the values oi x = (p^o and y = jlo — 1 which yield the minimum binding energy 
can be approximately solved. These values, we denote them as Xm = <^*om and ym = /tom — 1, 
is given by 

Xm = V*Om ~ (^24(7* (l - ^-P"*)) , Vm = hm " 1 ~ {?>C^f^^ . (63) 

These values lead to the mass Mj^, particle number of bosons NBm-i and particle number 
of fermions A^^m of the boson-fermion star of the minimum binding energy. The approximate 
values of them are 

1 C^" 

GMim ~ - + 



8 o^S/Tr/sr^ , /-, 3rn\^/'^^^ 



GmBNBm 



32^ [If, + (l - |F, 

(1 - 1^.)^'^ 



4 |F. + (l-|F. 



2/3 



3^ 

GmpNFr^ ^ * ^,^3. . (64) 



4 (|F. + (1 - 



One can sec that GmENsm + GmpNpm = 1/4 in this case. It is worth noting that the BTZ 
mass of the boson-fermion star approaches zero in the limit of C* = 0. 



VII. DISCUSSION 



We have obtained exact solutions describing boson-fermion stars with very large self- 
coupling constant in (2 -|- 1) dimensions. There is a critical value for C* = (A/(87rGm^))C, 
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where C is the absolute value of the (negative) cosmological constant. For C* > C* crit{F*) 
the binding energy of any solution is positive, the boson-fermion star configuration is not 
energetically favorable. 

To study the (2 + 1) dimensional boson star, with a finite value of the self-coupling or 
direct couplings between bosons and fermions, we must perform numerical analyses. The 
present study of the exact solution is the ground for studying the properties of the boson- 
fermion stars in general cases. 

For a finite self-coupling case, the external metric cannot be described by the exact BTZ 
metric. We expect that the external metric will approach the black hole metric in the large 
distance asymptotically if the global attractive force due to the cosmological constant is 
sufficiently effective. On the other hand, another types of the external metric will also be 
expected, which exhibit different behaviors from the BTZ metric. The possibility is worth 
studying. 

In the present paper, we have discussed only static configurations. We are interested 
also in the process of formation of the boson-fermion stars. The formation of stars or black 
holes (or possible naked singularities) will be clarified by the investigation of time-dependent 
processes. The quantum mechanical processes, including even the tunneling process, may 
play an important role in the creation of small boson-fermion stars. This subject will be 
considered in a separate research in the future. 

The time-dependent solution describing radial pulsations in (2 -|- 1) dimensions can be 
studied first of all. It will be discussed in elsewhere. 
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FIG. 1. Contour lines of constant Nb (solid lines) and Np (broken lines) in the {x, y) plane, 
where x = <^*o and y = jlo — 1. The bold line shows the line of zero eigenvalue (see text). Contour 
lines of the binding energy is also shown by dotted lines. In these figures, is set to unity, and 
C* is (a) 1, (6) 5, and (c) 9. 
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FIG. 2. The value of minimum binding energy of a boson-fermion star as a function of and 
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